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Infinite families of multi-indexed orthogonal polynomials are discovered as the solutions of exactly solvable one- 
dimensional quantum mechanical systems. The simplest examples, the one-indexed orthogonal polynomials, are the 
infinite families of the exceptional Laguerre and Jacobi polynomials of type I and II constructed by the present authors. 
The totality of the integer indices of the new polynomials are finite and they correspond to the degrees of the ' virtual 
state wavefunctions' which are 'deleted' by the generalisation of Crum-Adler theorem. Each polynomial has another 
integer n which counts the nodes. 
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1 . Introduction 



We will report on the discovery of infinite families 
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S of multi-indexed orthogonal polynomials, which form the 
main parts of the eigenfunctions of exactly solvable one- 
dimensional quantum mechanical systems. This will add 
£NJ a huge number of new members to the inventory of ex- 
actly solvable quantum mechanics 0, [2} , which had seen 
2$ rapid increase due to the recent discovery of various kinds 
of infinite families of exceptional orthogonal polynomials 
The method for deriving the multi-indexed or- 
thogonal polynomials is simple generalisation of the well- 
known method of Crum Q and its modification by Krein- 
Adler Starting from an exactly solvable Hamiltonian 
t-H system, Crum's method provides another exactly solvable 
^ system by successively deleting the lowest lying eigen- 
• <-h state. The modification due to Krein-Adler uses dele- 
tion of a finite number of higher eigenstates satisfying 
^ 'certain conditions. The present method is based on the 
"same spirit of 'deleting' a finite number of virtual states 
instead of the eigenstates. Here the virtual state wave- 
functions are the solutions of the Schrodinger equation of 
the original system and they have negative energies, with 
respect to the groundstate eigenvalue. The virtual state 
wavefunctions have no node in the interior and they are 
square-non-integrable. Their reciprocals are square-non- 
integrable, too. We apply this method to the Hamiltonian 
systems of the radial oscillator potential and the Darboux- 
Poschl- Teller potential, which are the best known exam- 
ples of exactly solvable quantum mechanics [l|, 0] . Their 
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eigenfunctions consist of the Laguerre and the Jacobi poly- 
nomials, respectively. 

When only one virtual state of degree £ (> 1) is deleted, 
the resulting system is the exceptional Xg Laguerre or Ja- 
cobi polynomials of type I or II [H, [Bj] . In general M dis- 
tinct virtual states of degree Vj (> 1) are deleted and the 
resulting eigenfunctions consist of orthogonal polynomials 
indexed by the integers (v%, . . . ,v«). The multi-indexed 
orthogonal polynomials are natural generalisation of the 
exceptional orthogonal polynomials. They start at degree 
£ > 1 instead of a degree constant term, thus avoiding the 
constraints by Bochner's theorem [ll|, [lij]. As solutions of 
exactly solvable quantum mechanical systems they consti- 
tute a complete system of orthogonal functions. The first 
examples of the exceptional Laguerre and Jacobi polyno- 
mials are constructed by Gomez-Ullate et al. [l2[ in the 
framework of Sturm-Liouville theory and by Quesne [l3| 
within the shape invariant 14] and exactly solvable quan- 
tum mechanics. They are equivalent to the lowest mem- 
bers (£ = 1) of the Xg Laguerre and Jacobi polynomials 
derived in [3j by the present authors. The present method 
could be considered as the generalisation of the Darboux- 
Crum transformations [lj| applied for the derivation of the 
exceptional orthogonal polynomials [l6|-(20l]. 

It is expected that these new polynomials would play 
important roles in physics, mathematics and related dis- 
ciplines. The multi-indexed Jacobi polynomials provide 
infinitely many global solutions of second order Fuch- 
sian differential equations with 3 + £ regular singulari- 
ties 21] . Other theoretical developments are: non-linear 
identities underlying the shape invariance |22j | are clari- 
fied, the structure of the extra zeros of the new orthog- 
onal polynomials l'2-j are exemplified, two-step Darboux 
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transformations 24| are explicitly constructed for the type 
II exceptional Laguerre polynomials, various use of non- 
physical states 25l|26| for generating the exceptional poly- 
nomials are suggested and several applications 271 l28| are 
reported. In this connection, let us mention references 
of multi-step Darboux-Crum transformations, higher or- 
der susy transformations, A/"-fold susy transformations, etc 
[29]. Some of them utilise non-physical solutions. 

This Letter is organised as follows. The general scheme 
of virtual states deletion is presented in section two. The 
explicit forms of the new polynomials are derived in section 
three, together with the forward-backward shift operators 
and the second order differential operators. The final sec- 
tion is for a summary and comments. 

2. Adler-Crum Scheme for Virtual States Deletion 

Here we present the method of virtual states deletion 
in its generic and formal form. For rigorous treatments of 
possible singularities at the boundaries, we need to specify 
the details of the Hamiltonian systems, which will be done 
in the next section. Throughout this paper we choose all 
the solutions of the Schrodinger equations to be real. We 
start with an exactly solvable one dimensional quantum 
mechanical system defined in an interval x\ < x < x 2 , 



H(j> n (x) = £ n (j> n {x) (neZ> ), 



dcf 



d,X 4) m {x)* 4> n {x) = h n 5r, 



(1) 

(2) 



For simplicity, we assume that it has discrete eigenvalues 
only with vanishing groundstate energy: = £0 < £\ < 
£ 2 < • • • , which allows us to express the potential in terms 
of the groundstate wavefunction (f>o(x) which has no node: 



(3) 



Let us suppose that the Schrodinger equation has a finite 
or infinite number of virtual state solutions with negative 
energies: 

H4> v {x) = SvMx), £ v <0 (v e V), (4) 

which have no node in the interior x\ < x < x 2 . Here V 
denotes the index set of the virtual states. We also require 
the infinite norm conditions, (<p v ,(p v ) — (1/</> v , 1/0 V ) = 00 
and this means that the virtual state wavefunctions vanish 
at one boundary and diverge at the other: 



lim <p v (x) — lim d x <fi v (x) = 0, 



lim <p v (x) — lim d x 4> v (x) = 00, 

X^X^ X^Xfr 



(5) 



where j^fce{l,2} and Vv 6 V. The set of eigenfunctions 
{(f> n (x)} forms a complete set in an appropriate Hilbcrt 
space H, for example: 

H = f {4>(x) e L 2 [xi 1 x 2 ] I lim (a; - xi)^d x tp(x) = 0, 



lim (#2 — x)?d x i/j(x) = 0}. 

xfx 2 

(For xi — —00 and x 2 = 00, replace (x — x\)? — > (~x)i 
and (x 2 — x)? — > x? .) That is any function in H which is 
orthogonal to all the eigenfunctions has a zero norm: 



(tf,&0=0 (VneZ> ) 



= 0. 



After M-steps of deleting the lowest lying eigenstates 
by Crum's method the resulting Hamiltonian system 
is iso-spectral to the original one given by 

(x) = £ n 4™ ] (x) (n = M, M + 1, . . .), 

W[0 O ,0l, ■ • . ,(f>M-l]{x) 
M-l 

3=0 

JjlM] {x) dcf u{x) _ 2 q2 log ]( W[0O) ^ . . . ; ^ M _ l]{x ) I , 

in which W[/ x ,--- ,/„](*) d =? f det(^- 1 /fe(a:)) 1 < iifc < n isthe 
Wronskian. These formulas have been repeatedly rediscov- 
ered under various names, higher derivative susy, TV-fold 
susy, etc j29[. The corresponding Adler-Crum formulas 
for deleting the set of M distinct eigenstates specified by 
V = {d\,d 2 , ■ ■ ■ ,d M }, dj > 0, are 

rt M U [ n 1] (x) = £ n ^ ] (x) (n e Z> \27), 



,[M}( \ dcf W[0 dl , (j)d 2 , ■ ■ ■ , 4>d M , <Pn] (x) 



W[</» < j I ,0d 3 ,...,^ Ar ](a:) 



M 



(^\^) = Y[(£ n ^£ d] )-h n S mni 

3=1 

TJ[M] {x) d|f u{x) _ 2Q 2 l 0g \ W[(f)di , d2 , . . . , faj ( X ) 



The set V must satisfy the conditions rX/=i( TO — ^3) — 
0, Vm £ Z>o [13, [33]. The exactly solvable Hamiltonian 
system obtained by deleting M distinct virtual states T> = 
{d\, . . . , djv/}, dj S V looks almost the same: 



«^( I )=£ n ^( I ) (n£Z> ), 
H^^(x)=£ v ^(x) (v£V\5), 

W[^ rfl ,0d 2 ,---,0d M ](^) 

M 

(<f>W \ = l[(£ n - £ dj ) ■ h n S mn , 
3=1 

l[M] ( \ dcf W[0 dl ,0d 2 ,...,^d M ,0v](a:) 



<t>r l (x) 



W[^d 1 ,^ a ,...,^ M ](a:) 



(6) 
(7) 

(8) 
(9) 
(10) 



JJ[M] {x) d£* jj {x) _ 2 q2 log | W [^ rfi , d2 , . . . , $ dM ]( x ) I . (11) 

It should be stressed that the final results after M deletions 
are independent of the orders of deletions. The algebraic 



side of the derivation of (j6j)- (TTTj) is basically the same as 
that for the Adler-Crum case. The essential difference is 
the verification of the inheritance at each step of the prop- 
erties of the virtual states; no-nodeness and the proper 
boundary behaviours. This is guaranteed by choosing the 
parameters of the theory g and h large enough, see ([23| 
and (|24|) . 

Here we show the process % — > "H' 1 '. It is elementary to 
show 



H = A\A d + £ d , 



2 def d 

Ai = -y 

ax 



d x (j)d{x) 
(j>d{x) 



yd e V. (12) 



Since (f> d (x) has n0 node, the operator Ad and A d are reg- 
ular in the interior. Next we define T-C 1 ^ by interchanging 
the order of these two operators: 

n^=A d Al+£ d , U^(x) = U(x)~2d 2 x \og\Mx)\. (13) 

It is straightforward to show that H and T-L^ are iso- 
spectral 



H^(x) = £ n ^(x), 
^\x)^AdMx)^ dAn]ix) 

n^(x) = s v ^(x), 



<t>d{x) 



(neZ 



4> [1] {x) d = AdMx) 



VJ[4> d ,4> v ](x) 



4>d{x) 



(v € V\{d}), 



4>n{x) 



A 



f 

L n 



Od 



4$) = (£n - £d)((f>m,(f>n) = {£n ~ £d)h n 8 mn . 

Obviously the state d is deleted from the virtual state spec- 
trum of %M , for Ad4>d{x) — 0. Since the zero-modes of A d 
and ,4^ are (f> d (x) and l/ck> d (x), which have infinite norms, 
the sets of eigenfunctions {(j> n {x)} and (a;)} are iso- 
spectral and in one to one correspondence. Thus the new 
set {0n'(2:)} is complete. If not, there must exist a finite 
norm vector \P in the Hilbert space which is orthogonal to 
all the basis: (*, ^L 1 ') = 0, Vn e Z> . This would lead to 
a contradiction {Al^, (j> n ) = 0, V?i G Z> , since A a ^ ^ 0. 
Next we show the logic for demonstrating that the virtual 
state solutions {(f>v\x)} have no node in the interior. By 
using the Schrodinger equations for them we obtain 

d x W[4> d ,fa](x) = (£ d - £ v )Mx)Mx), (14) 

which has no node. Since we can verify that W[(f> d , (j)v]{x) 
vanishes at one boundary, no-nodeness of W[4>d, 4>A(x) 
in the interior follows. When the virtual states d and 
v belong to different types (I and II) the Wronskians 
might not vanish at both boundaries. In that case we 
can show that they have the same sign at both bound- 
aries W[(j>d, <j>v](xi)W[(i)d, <j>v](x2) > 0. Likewise we can 
show, for all the explicit examples in the next section that 
<^>v' and l/$v have infinite norms. The steps going from 
— > and further are essentially the same. 



3. Multi-indexed Laguerre and Jacobi Polynomials 

Here we will apply the method of virtual states deletion 
to two well-known systems of exactly solvable quantum 
mechanics with the radial oscillator and Darboux-Poschl- 
Teller potentials: 



U(x) 



5(5-1) 



(l + 2.g), 



xi = 0, x 2 = oo, g> \, 
g(g-l) , h(h-l) 



(15) 



• 2 

sm x 



— (g + h) 2 , 

cos^ x 

x\ = 0, x 2 = § , g, h > |, 



in which L and J stand for the names of their eigenfunc- 
tions, the Laguerre and Jacobi polynomials. For the ac- 
tual parameter ranges consistent with specific deletions, 
see (|2"3")) - (f2"4")) . The eigenfunctions are factorised into the 
groundstate eigenfunction and the polynomial in the sinu- 
soidal coordinate r\ = T)(x) |3lj ]: 



(j> n {x; A) = <p (x; \)P n (r](x); A), 



(16) 



in which A stands for the parameters, g for L and {g, h) 
for J. Their explicit forms are 



P n { m g)^L { r h \nl 



I \ def 2 
TJ[X) = X , 



£ n {g) = 4n, 

def 



K{g) d = ^r r ( n 



I), 



clef 



(17) 



4>o(x;g,h) = (sin x) ff (cos x) rj{x) — cos2a; 



P n { m g,h) = Pi B ~ 



fa), 



dof 



£n(g, h) — An(n + g + h), 

h n(g,h) = — ■ - ■ ■ - ■ , v (18) 



2n!(2n- 



h)T{n 



They have two types of virtual states, which are again 
polynomial solutions. The total number is finite depending 
on the parameter values g and h, except for LI which has 
infinite. The virtual states wavefunctions for L are: 

LI: &(z)^e**V£fa(s);s), 

£l(v;g) = Pv(-rr,9), veZ> , 

^ f -4(. 9 + v+I), ~S Id ^-l, (19) 

L2: ^(^c^it'tloWij), 



o,i,..., [g -§]', 



(20) 



in which [a]' denotes the greatest integer less than a and 

5 ' will be used later. For no-nodeness of £*' n , see (2.39) 
of 22] . The virtual states wavefunctions for J are: 

Jl : $,{x) A = {sin x) 9 {coax) 1 " 11 i\{r,{x)]g,h), 



Clivi 9,h) = P v (r);g,l-h), v = 0,l,. 



.,[h- 



-1' 

21 ' 



I). 



5 Id = (-1,1), (21) 



4 1?-4(9 + v+ \){h- 
J2 : ^(z) d =S f (sin^^cosaO^ 1 ^);^), 

&{mg, h ) d = P v (m 1 ~9,h), v = o,i, . . . , [# 

-4( 9 -v-±)(M 



2J ' 



£ 11 = f 



5 IId =i f (l,-l).(22) 



The larger the parameter g and/or h become, the more 
the virtual states are 'created'. The LI system is ob- 
tained from the Jl in the limit h — » oo [5J. This ex- 
plains the infinitely many virtual states of LI. Let us 
denote by V I,n the index sets of the virtual states of 
type I and II for L and J. Due to the parity property 



of the Jacobi polynomial Pi a ' 0) (-x) = (-l)"P r ^'" ; (x), 
the two virtual state polynomials and £~ for J are re- 
lated by £v ! ( — Vi 5> h) = (— l) v Cv( r ?i h, gY For no-nodeness 
of $ n , see (2.40) of [H and (3.2) of @._It is also obvi- 
ous that the virtual state wavefunctions cf> v and their re- 
ciprocals 1/4> V are not square integrable in all four types. 
There is a third type of negative energy polynomial so- 
lutions for both L and J. For J, the degree n solution 

is (sinx) 1 ~ £ '(cosa;) 1_ ' l i^/ l 9+2 ' h+2 \n). Even when they 
have no node in the interior, their reciprocals are square 
integrable and they cannot be used for the virtual states 
deletion. The label is special in that the wavefunctions 



>(/3,«)/ 



satisfy <ft {x\ X)<Pq(x\ X) = c 



-1 dr)(x) 
dx 



since £j' n = 1. Here 
-4 for J. We will not 



the constant = 2 for L and c T 
use the label states for deletion. 



A 



■di 



A 



A. 



3 
2 
1 


hi 

d 2 " 
: "i 



H 

original d\ deleted 
system system 



A\ 



d x di 



y 



d\ , d 2 deleted 
system 



M+N virtual states 
deleted system 



Figure 1: Schematic picture of virtual states deletion. 
The black circles denote eigenstates. The down and up 
triangles denote virtual states of type I and II. The 
deleted virtual states are denoted by white triangles. 

Since there are two types of virtual states available, the 
general deletion is specified by the set of M + TV positive 



,d l Ml d l ^...,d 1 ^}, d j > which 



integers V = {d\, . . 
are the degrees of the deleted virtual state wavefunctions. 
Note that the subcases of either M = or N = are 



meaningful. In order to accommodate all these virtual 
states, the parameters g and h must be larger than certain 
bounds: 



L : g > max{iV + §, d" + |}, 
J: g > m&x{N + 2,df + ±}, 
h > max{M + 2,d) + |}. 



(23) 



(24) 



See Figure 1 for the schematic structure of virtual states 
deletion. 

Like the original eigenfunctions (fTBf the n-th eigenfunc- 
tion 4>x>, n {x;\) = (ft^'^ix) after the deletion d9]) can be 
clearly factorised into an x-dependent part, the common 
denominator polynomial Sx> in r\ and the multi-indexed 
'polynomial Px>,n i n V- 



4>^ N \x) = ^ <n (x;X) 
ipT>(x;X) 



c^ [+N ^{x-X)P v ^{r 1 {x)-X), 



def (f>o{x;X [M:N] ) 
E v (rj(x);X) ' 



(25) 



in which (J)q(x;X) is the groundstate wavefunction (jTTjl - 
(fT8"|), Here the shifted parameters X [M ' N] after the [M,N] 
deletion are X [M ' N] d = A - MS* - NS™, explicitly they are 



X [M,N] = 



g + M - N 
{g + M -N,h 



M + N) 



L 
J 



(26) 



Needless to say, the denominator polynomial has no 
node in the interior. The polynomials Px>, n and are 
expressed in terms of Wronskians of the variable r\: 



Pv,n(W>X) = W[/ii, . . . ,Hm,Ui, . ..,VN,Pn](v) 



e -Mt) „(M+g+|)iV 



rL-r,\ (M+g+^N (l±n_\ (N+h+±)M 



Bp^AJ^WtMi,.. 



.,Hm,vi,---,^n](v) 



e -Mr) „(M+ fl -|)AT 



l-r,^ (M+g-f )iV / i +v ^ (N+h-±)M 

3 



L 
J 

L 
J 
L 

J 

: L 
: J 



(27) 



(28) 



(29) 



in which P n in (|27p denotes the original polynomial, 
P n (v'i 9) f° r L and P n (?7; g, ft.) for J. These polynomials de- 
pend on the order of deletions only through the sign of 
permutation. The multi-indexed polynomial Pv.n is of de- 
gree £ + n and the denominator polynomial S-p is of degree 
£ in 7?, in which £ is given by 



M N 

1 = ^4+^4 I -±M(M-l)-~/V(A/'-l)+MiV > 1. 



(30) 



4 



Here the label n specifies the energy eigenvalue £ n of 4>T>,n- 
Hence it also counts the nodes due to the oscillation the- 
orem. Although they do not satisfy the three term re- 
currence relations, the multi-indexed polynomials {Px>,n} 
form a complete set of orthogonal polynomials with the 
orthogonality relations: 



(X)6 

ii in 

Ylf=xin + 9 + d} + l)-UU(n + 9-df-^ :L 
< A- M ~ N njii(n + 9 + d) + |)(n + h-d)-\) , 
x nf=i(« + 9 - df - ±)(n + h + df + |) : J 

(31) 

where the weight function of the original polynomials 
A)c?r/ = (j)o(x; \) 2 dx reads explicitly 



Wfa A) 



clef 



23 + 



^(1 -77)3-3(1 + r)) h ~2 



L 

J 



(32) 



This is obtained by rewriting the orthogonality relations 
of the eigenfunctions ([§]) after the [Af , N] deletion. 

In the rest of this section we explore various proper- 
ties of the new multi-indexed polynomials \Py.n\- As for 
the exceptional orthogonal polynomials 0, 0, H, EH 0, 0] j 
the lowest degree polynomial Px> i o(^;A) is related to the 
denominator polynomial Sd(?7; A) by the parameter shift 
A -> A + 8 (8 = 1 for L and 8 = (1, 1) for J): 



Pt>,o(v,X) =Sx,(t7;A + 8) 



(-i) M nf=i(<? 



d 11 



I) 



x(-2)- w nL(3 



4 



(33) 



d 11 



The virtual state wavefunction with the index v G T> lead- 
ing to the final [Af , N] deletion is 4>[ M ' ' N ' ] with the 
index set V d = V\{v}, in which [Af, N'] = [Af - 1, N] if 
v G V 1 , and [Af, N - 1] if v G V 11 . It has the following 
form: 



I[M',JV' 



Or; A) 



jtf+iy-i Sp(r?(a;); A) 



x < 



(-i) N M^\ lM ' 1 ' N] )i- {9+M - N - 1 \ v e v 1 

^(^^(A^^- 1 !)), v € V 11 

(-l^oM 1 ^- 1 ^)), veV 

0o(&; t n (A' M,JV -^)), vev" 



L 

L 
J 
J 

(34) 



Here the twist operator t acting on the parameters is de- 
fined as t u (g) = 1 — g for L, t l (g,h) = (g, 1 — h) and 
t u (g, h) = (l-g, h) for J. Note that the virtual states ([20])- 
(22) have the form <j> v (x; A) = </> v (a;; t(A)). These relations 
are essential for the determination of the lower bounds of 



the parameters given in (f23l) - (j2"4|) . which guarantee that 
the virtual state properties are correctly inherited at each 
step of deletion. 

The Hamiltonian H v (A) = H [M ' N] of the [Af, N] deleted 
system can be expressed in terms of its groundstate eigen- 
function <jyp o(x; A) = 0Q M,Af '(a;; A) with the help of (l33|) : 

. def d d x <j)T>.o(x; A) 
Av{X) = -j t — t — tt-, 

dx (pT>,o{x;X) 

a i w a ( ^\M.N] ^v(v(x); A + 8) 
(f>v,o( x i A) oc (po{x;X l ') - . . , . 

£z v (r)(x); X) 



(35) 



(36) 



This has the same form as the various Hamiltonians of 
the exceptional orthogonal polynomials including those of 
the discrete quantum mechanics d, 3, [E Bl ■ Reflecting the 
construction 0, @] it is shape invariant [3, S [H| 

Av(X)A v (X) f =A v (X + S)^Av(X + 8) + £ 1 (X). (37) 

This means that the operators At>(X) and AviX)^ relate 
the eigenfunctions of neighbouring degrees and parame- 
ters: 

A v {X)(f>vA^ A) = /n(A)0D, n -i(x; A + 8), (38) 
Av(Xy<t)v,n-i(x; X + 8) = 6„_i(A)0 23 ,„(a;; A), (39) 

in which the constants / n (A) and 6 n _i(A) are the factors 
of the eigenvalue / n (A)b n _i(A) = £„(A): 

/„(A) = { Zl {n + g + h) !j. 6»-i(A) = -2n :L,J. 
The forward and backward shift operators are defined by 
T V {X) d = *Pv (x; X + 8)- 1 o A V {X) o ^ p (a; A) (40) 

C " Sp(r7;A) Ur? -2,(77; A + 5) /' 1 ' 
B V (X) d = i>v {x; Xy 1 o Av{X)^ o ^ (x; A + <5) (42) 



-4c^ 1 c 2 (r?)- 



S c (t7;A + ^) 



C2(r?) 



Sp(?7; A) 



, (43) 



in which the functions ci(f];X) and C2(r)) are those ap- 
pearing in the (confluent) hypergeometric equations for 
the Laguerre and Jacobi polynomials 

L: a(r), A) = f g + \ - 77, 02(77)=%, 

J: ci(t7, A) d =h-g- (g + h+l)r), c 2 (77) = f 1 - t? 2 . 

Their action on the multi- indexed polynomials Pv,n{v\ X) 
is 

Fv{X)P v ,n{mX) = f n (X)P v , n - 1 (r ] ;X + 8), (44) 
B v (X)Pv,n-i(v,X + 8) = b n ^(X)P v , n (r);X). (45) 



The second order differential operator "Hp (A) governing 
the multi-indexed polynomials is: 

n v (\) d = ipr>(x; A) -1 o %v (A) o ^{x- A) = B V {\)F V {\) 

._ 4 ( C2( , ) ^ +(ci( „ A ^ ) _ 2c2( , ) ^) ) j. 

d%Z v (r);\) x [m,n] cn^S-d^; A) 

+ C2W)-h- — r- ci(r/, A L 1 - d)-^- 



Bp(?7; A) 



Sp(??; A) 



(46) 
(47) 



Since all the zeros of £15(77; A) are simple, ([4"Tf is a Fuchsian 
differential equation for the J case. The characteristic ex- 
ponents at the zeros of St>{tj; A) are the same everywhere, 
and 3. The multi-indexed polynomials {Pjy t n(r); A)} pro- 
vide infinitely many global solutions of the above Fuchsian 
equation (|4"Tj) with 3 + £ regular singularities for the J case 
[2lJ. The L case is obtained in a confluent limit. These sit- 
uations are basically the same as those of the exceptional 
polynomials. 

Although we have restricted d*' 11 > 1, there is no ob- 
struction for deletion of d 1 ^ 11 = 0. In terms of the multi- 
indexed polynomial (|^T)) . the level deletions imply the 
following: 



Pv>, n (m*-s ) x A 



1,4* + !,...,< + !}, (48) 



ffe'.nfaA-a ) X B 



1. 



d u 

■ ! a N-l 



Pv,n(ji;\) 

V = {d\-l,...,d I M _ 1 

V = {d\ + l,...,d 1 M 
where the multiplicative factors A and B are 

A=l -(-2)- M nf = T 1 1 (9-h + d) + l) 

^(-^- N Uf=M 1 + l)-(n + h-l) 



1}, (49) 



B 



(-l) M U/=M) + l)-(n + g-l) 

2- M n-Ii(4 + 1 )-(- 2 ) 



J 
L 



— N 



l 3 

x n 7 =i (k 



+ 4 I + l)-(n + . 9 -i) 



From (|33p . Sp behaves similarly. Therefore including the 
level deletion corresponds to M + N — 1 virtual states 



deletions. This is why we have restricted d 1 - 11 > 1. 

These relations (|4"5 |) - (|4l?|) can be used for studying the 
equivalence of 'Hp. It should be stressed that the same 
polynomials can have different but equivalent sets of multi- 
indices 2?, which defines the denominator polynomial Sp. 
Since the overall scale of Sp is irrelevant (see (131))) ') for the 
multi-indexed polynomial {-Pp, n }, we consider the equiv- 
alence class of Sp. For example, it is straightforward to 
verify the following equivalence and their duals I o II: 

S {1 i,..., fc i } (7 7 ;A-$ n ) -H^njfoA-fcS 1 ) (fc > 1), (50) 



J {m I ,...,(fc+m) 1 } 



n (77; A — md 



~ -{(fc+i)",...,(fc+m) n }( ? 7; A - ( k + l ) s ) ( fc > !)• ( 51 ) 

The simplest, fc = 1, of (15U1) corresponds to the well known 
fact that the X\ type I and II polynomials are identical for 
L and J Deleting the first k consecutive excited states, 
like the l.h.s. of (fST))) has led to many interesting phe- 
nomena in the ordinary and discrete quantum mechanics 
30|, 33, 32 1 . Classification of the equivalent classes leading 
to the same polynomials is a challenging future problem. 

The exceptional X p o rthogonal polynomials of type I 
and II, [H [H, i, i, M [HI correspond to the simplest 
cases of one virtual state deletion of that type, V — {i 1 } 
or {£ n }, I > 1: 

tt( m \)=3 D (ri;\ + £6 + S), (52) 
P i , n (r } ;X) = P- D , n (r];X + £S + S) x A, (53) 



~i,ii 



and the multiplicative factor A is A = — 1 



lw for XL2, 2(n + ft + A)" 1 for XJ1 



2 / 
f)" X 



for XJ2. Most formulas between 



where S = S 
for XL1, (n + g 
and —2(n + g -\ ^ 
([55]) and (|4"Tj) look almost the same as those appearing 
in the theory of the exceptional orthogonal polynomials 

S I, S El [3 0,10. 

In a recent paper [24[ Gomez-Ullate et al discussed "two- 
step Darboux transformations," which corresponds to the 
examples of T> = {m^mi, 1 } for L in our scheme. This 
paper stimulated the present work. 

4. Summary and Comments 

This is a first short report on the discovery of infinite 
families of multi-indexed orthogonal polynomials. They 
are obtained as the main part of exactly solvable quan- 
tum mechanical systems, which are deformations of the 
radial oscillator and the Darboux-Poschl- Teller potentials. 
Although these new polynomials start at degree I (> 1), 
they form a complete set of orthogonal polynomials. The 
exactly solvable modified Hamiltonian and the eigenfunc- 
tions are obtained by applying the modified Crum's theo- 
rem due to Adler to the original system to delete a finitely 
many virtual state solutions of type I and II. The sim- 
plest case of one virtual state deletion reproduces the ex- 
ceptional orthogonal polynomials [HI, O d, HI [Hj]. For 
the Jacobi case, the new polynomials provide infinitely 
many examples of global solutions of second order Fuch- 
sian differential equations with 3 + 1 regular singularities 
2l| . Like the undeformed theory, the Laguerre results can 
be obtained from the Jacobi results in certain confluence 
limits. But we presented all results in parallel, for better 
understanding of the structure. The parameter ranges (|23p 
and (|2~4"]) are conservative sufficient conditions with which 
the proper boundary behaviours and the completeness are 
guaranteed. 

As a first report, only the basic results are presented. 
Some important issues cannot be included due to space 



restrictions, for example, an alternative proof of shape in- 
variance, the action of various intertwining operators, etc. 
Many important problems are remaining to be clarified. 
For example, the bispectrality [l^|, generating functions 
21fl , properties of the zeros of the denominator polynomial 
Sx> and the extra zeros of the multi-indexed polynomials 
{Pv,n} HH, classification of the equivalent classes leading 
to the same new polynomials. 

The exceptional Wilson, Askey- Wilson, Racah and q- 
Racah polynomials were constructed in the framework of 
discrete quantum mechanics 0, [2(| 0, [1] , and the Crum- 
Adler theorem for discrete quantum mechanics were also 
presented [33, 

Hoi S3. 

The method of virtual states dele- 
tion presented in this Letter is applicable to discrete quan- 
tum mechanics and it is easy to write down general for- 
mulas like as (joT)- (fTTj) . This gives multi- indexed orthogo- 
nal polynomials corresponding to the deformations of the 
Wilson, Askey- Wilson, Racah and g-Racah polynomials. 
In concrete examples, more than two types of virtual state 
solutions are available and much richer structures are ex- 
pected. We will report on these topics elsewhere. 

Before closing this Letter, let us make two small com- 
ments. The original Schrodinger equation ([lj has a general 
solution for an arbitrary energy E, T-Ltp{x) — £ip(x). If we 
define the corresponding [M, N] deletion solution 

- [M,N] f„\ def W[<^ , . . . , , ftl ,...,4> I J N ,^} (x) 



W[4> J dl ,...,4> l dM ,4>l,...,4> l l](x) ' 



it solves the deformed Schrodinger equation with the same 
energy, too 

n [M.N]^lALN] {x)=£i; [M.N] {x y 

Among the multi-indices, we can include a finite number 
of continuous real parameters as labe ling the type I vir- 
tual solutions as described in 0, llVj The resulting 
orthogonal functions are no longer polynomials, but most 
properties of the multi-indexed polynomials are shared by 
these new orthogonal functions. 
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